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Abstract 


A mathematical model was developed for modelling the performance of solid oxide fuel cell (SOFC) with functionally graded electrodes at the 
micro-scale level. The model considered all forms of overpotentials and was able to capture the coupled electrochemical reactions and mass transfer 
involved in the SOFC operation. The model was validated by comparing the simulation results with experimental data from the literature. Additional 
modelling analyses were conducted to gain better understanding of the SOFC working mechanisms at the micro-scale level and to quantify the 
performance of micro-structurally graded SOFC. It was found that micro-structural grading could significantly enhance the gas transport but had 
negligible effects on the ohmic and activation overpotentials, especially for thick electrodes. However, for thin electrodes with large particles, 
too much grading should be avoided as the increased activation overpotentials may result in higher overall overpotentials at a medium or low 
current density. Among all the cases tested in the present study, the micro-structurally graded SOFC showed significantly higher power density than 
conventional SOFC of uniform porosity and particle size. The difference between micro-structurally graded SOFC and conventional SOFC is more 
pronounced for smaller electrode—electrolyte (EE) interfacial particles. Particle size grading is generally more effective than porosity grading and it 
can increase the maximum power density by one-fold in comparison with conventional SOFC. The present study reveals the working mechanisms 


of SOFC at the micro-scale level and demonstrates the promise of the use of micro-structural grading to enhance the SOFC performance. 


© 2007 Elsevier B.V. All rights reserved. 
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1. Introduction 


Solid oxide fuel cells (SOFCs), which are efficient, tolerant 
of impurities, and flexible in fuel choices, will play an impor- 
tant role in the next generation of energy technology [1-3]. 
In order to further enhance its energy conversion efficiency, 
efforts have been made in developing new materials and design- 
ing novel geometries. Hwang et al. [4] developed a numerical 
model to investigate the characteristics of a mono-block-layer- 
built (MOLB) SOFC and they found that the MOLB-SOFC had 
a higher fuel/oxidant utilization than the planar SOFC. Ramakr- 
ishna et al. [5] performed numerical analyses and proposed an 
innovative thin-walled geometry to improve the power density of 
SOFC. More recently, a novel 3-D SOFC was developed by Koh 
et al. [6] to increase the volumetric power density by increasing 
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the volumetric surface areas of the electrodes. The philosophy 
of these studies is to optimise the fuel/reactant flow field at the 
macro-scale level so as to improve the cell performance. 
Functionally graded electrodes (FGE) have been applied to 
SOFC in the recent years and offered potentials to enhance the 
SOFC performance at the micro-scale level [7—10]. Grading can 
be classified into three ways: (1) composition grading, (2) parti- 
cle size grading, and (3) porosity grading. The porosity grading 
and particle size grading are two forms of micro-structural 
grading. The aim of grading is to optimise the electronic/ionic 
conductivity, to increase the electrochemical reactivity of the 
electrode, or to enhance the gas transport. Despite some exper- 
imental works available in the literature, the modelling studies 
on FGE are limited. Schneider et al. [10] investigated the per- 
colation effects in composition graded SOFC electrodes with 
Kirchhoff’s current law. It was found that the graded elec- 
trodes with varying composition did not perform significantly 
better than the non-graded electrode but the percolation was 
improved. Greene et al. [11] developed a macro-scale numerical 
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model to investigate the mass transport phenomena of poros- 
ity graded SOFC electrodes. The effects of porosity change 
on the concentration and ohmic overpotentials were consid- 
ered in their macro-scale model but the activation overpotential 
was treated as a property independent of the porosity as the 
electrochemical reactions were assumed to occur merely at 
the electrode—electrolyte (EE) interface. With a uniform pore 
size along the depth, an enhanced SOFC performance (8.7% 
improvement) was observed for the porosity-graded electrodes 
due to the reduced transport resistance [11]. 

Despite the abovementioned investigations, the studies on 
micro-structurally graded SOFC are insufficient and its mech- 
anisms are not well understood. First, the literature is lacking 
detailed modelling of the particle size graded SOFC. In addition, 
the existing study on porosity graded SOFC at the macro-scale 
level is insufficient to provide details to explain how porosity 
grading affects the SOFC performance. The present study aims 
to analyse the performance of micro-structurally graded SOFC 
by micro-scale modelling. The porosity grading and particle 
size grading have been investigated and their potentials for use 
in SOFC have been evaluated by comparison with non-graded 
SOFC. The model developed is able to predict the micro-level 
distributions of gas concentration, current densities, and elec- 
trochemical reactions along the electrode depth as well as the 
macro-level SOFC performance. The results presented in this 
paper are valuable for design optimisation of SOFC. 


2. Model development 
2.1. Working mechanisms of SOFC 


The coupled electrochemical reaction and mass transfer in a 
non-graded composite anode are illustrated in Fig. 1. The reac- 
tion processes can be summarized as: (1) transport of reactant 
(H2) to the reaction sites through the pores of the anode and trans- 
port of O?~ from the electrolyte to the reaction sites through the 
ionic conducting particles; (2) electrochemical reaction of H2 
and O°- to form H20 and electrons at the active sites; (3) trans- 
port of electrons from the active sites to the current collector 
through the electronic conducting particles and transport of the 
H20 product to the anode surface via the pores of the anode. 
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(O electronic conducting particles 
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Fig. 1. Schematic of the composite anode of a SOFC, consisting of both elec- 
tronic conducting particles and ionic conducting particles. 
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Fig. 2. Schematic of the functionally graded anode: (a) porosity grading with 
constant particle size; (b) particle size grading with constant porosity. 


Compared with a conventional non-graded electrode (Fig. 1), a 
graded electrode can be formed by either one of the two differ- 
ent types of micro-structural grading, namely porosity grading 
and particle size grading. Fig. 2 shows the two types of grading 
done on an anode. Both porosity and particle size can vary along 
the depth of the anode, from a larger value at the anode surface 
to a smaller value at the EE interface. A graded cathode can be 
formed in similar manners. 


2.2. The anode 


Charge balance in the electronic/ionic conducting particles 
can be written as 


V. Je,a = —SvJna =-V. Jia (1) 


where Je a and Ji a is the current density (A m7?) in the electronic 
and ionic conductors, respectively; Sy the electrochemically 
active internal surface area per unit volume of the porous anode 
(m? m7), which depends on the electrode depth in an FGE; 
and Jn, is the transfer current per unit area of reactive surface 
(Am7?). 
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The active surface area (Sy ) can be calculated with the method 
based on the binary random packing theories [12-15] 


ZeZi 
Sy = msin?Or2nynen; i PoP; (2) 


where 0 is the contact angle between the electronic and ionic 
conductors; re the radius of the electronic conductors; m the 
total number of particles per unit volume; ne and nj; the number 
fractions of electronic and ionic conductors, respectively; Ze and 
Z; the coordination numbers of electronic and ionic conductors, 
respectively; Pe the probability of electronic conductors con- 
nected with the porous media; and P; is the probability of ionic 
conductors connected with the porous media. 

The total number of particles per unit volume, depending on 
porosity £ and particle sizes (radii re and r; for electronic con- 
ductors and ionic conductors, respectively), can be determined 
according to [14], 


l-e 
~ (4/3)ar3[ne + (1 — ne)(ri/ re) 


Other parameters used in Eq. (2) can be calculated by the same 
way with Costamagna et al. [13] and Chan and Xia [12]. In the 
present study, re is assumed equal to r; at the same electrode 
depth. 

The transfer current density Jn,a can be determined by the 
generalized Butler Volmer equation [13], 


tef Pu, az Fna 
na = Y0.a po exp RT 


Py (0) d = a)z Fn; 
g exp ( = (4) 
Pipo 


where he is the reference exchange current density; 
Pp,o and POs the partial pressures of H2O within the porous 
anode and at the anode surface, respectively; similarly, 
Px, and Pi the partial pressures of H2; œ the charge transfer 
coefficient, which is 0.5 for SOFC [16]; z the number of charges 
involved per reaction; and na is the overpotential in anode. The 
value of ņa is defined as [12,13] 


na = (G2, — boa) — (be, — Pia) (5) 


where ġia and @e a denote the ionic and electronic potentials, 
respectively; the superscript 0 represents the equilibrium condi- 
tion. 

The electronic and ionic potentials can be obtained by apply- 
ing Ohm’s law, 


(3) 


nt 


V - bea = p Jea and V. Qia = pile (6) 


where p°f is the effective resistivity (Q m); J the current density; 
and the subscripts i and e represents, respectively, the ionic and 
electronic conductors. The effective resistivity of electronic and 
ionic conductors can be determined by [14] 

E E 


eff eff 
——— d a — 7 
a= aeaa ™ ee a 


where @ is the volume fraction of the electronic conducting 
particles; £ the tortuosity of the electrode; cea and oj, are the 
conductivities of solid electronic conductor and ionic conduc- 
tors, respectively. 

The first and second derivatives of ņa can be written as 


dna ff ff 
dx = Paia n Pe aJe,a (8) 
and 
dn dJi a dJe, d 
Gi = Pla qe 7 Pea ge = Shala + Pera) 
P azF 
fr eff ff H Ma 
= Sv Jo a(i a T Pea [ae SRP ( RT ) 
l — F 
i exp ( aE: y) (9) 
Po RT 


At the anode surface, Jia=0 and Jea=Jtotal. At the EE 
interface, Jia =Jtota and Jeg =0. Therefore, the two boundary 
conditions for the above second order differential equation (Eq. 
(9)) can be written as 


dn, d 
x=0, = = — pe Jtotal and x= da, a 
dx x=0 x x=dy 
= pF total (10) 


In the steady state, the conservation of gas species can be 
written as 


V- Ni = Ri (11) 


where N; denotes the rate of mass transport (mol m~? s7!) 
and Ri represents the rate of reaction inside the porous elec- 
trode (mol m~? s~!). The rates of reaction for H2O and H3 can 
be expressed as Ry, = — Sy Jn,a/2F and Ry,0 = SvJn,a/2F, 
respectively. 

The mass transfer process of Hz in the porous electrode is 
by means of diffusion. Considering both Knudsen diffusion and 
molecular diffusion, the mass flux of hydrogen can be deter- 
mined by the Dusty Gas Model (DGM) as 


a (12) 


-—1 
1 — ayn, dyn, 
ff 
Hy, k Dio- 


P 1 
Np, = RT pet 


where æ = 1 — ,/My,/Mu,0; Mp, and Mp,o the molecular 
weights of Hz and H20, respectively; and De and DË jii 
are the Knudsen diffusion coefficient of Hz and the binary 
diffusion coefficient of H2—-H20 mixture, respectively. The cal- 
culation of the abovementioned diffusion coefficients, which are 
functions of the microstructures of the electrodes, can be found 
in references [17,18]. The important pore radius can be obtained 


by the expression developed by Nam and Jeon [19] as 


2 TeE 
Tp => 
P™ 31—e 


(13) 
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Substituting Eq. (12) into Eq. (11), the controlling equation 
for H? diffusion becomes 


dym, P a 1 1 — ayn, (Se ) 
2 ff ff ff 
dx Dii,0-H Diy k Di O-H; dx 
Sy Jaa RT 1 1- 
AFP eff go = (14) 
Dik Di,0-H 
The boundary conditions for Eq. (14) are given below 
mle 0 . dyp, = RT Jiotal 1 1 — ayn, 
= yn, = 
á > dx fzo 2FP \ Di, Dom 
(15) 


In the abovementioned equations, the gas molar fraction, gas 
diffusion coefficients, electronic/ionic current density, reaction 
rate, and the overpotential are dependent of depth x. Due to 
the complicated phenomena involved, an iterative scheme is 
used to solve the non-linear equations. The iterative compu- 
tation is successful if a convergence is obtained. After obtaining 
the distributions of species concentration, current densities, and 
overpotentials, the overall electrode overpotential can be deter- 
mined [13], 


Nea = (Ola — Poa) — Pe,a(x=0) — Pisa(x=da)) (16) 
2.3. The cathode 
Similar to the anode, the coupled electrochemical reactions 


and transport phenomena at the SOFC cathode can be deter- 
mined, 


= Sy JpeCofe + pee) 7 exp (S) 
a ( a =o) | (17) 
and 
V- Jee = —SvIne = —Sv IS py i (=) 
ah ( a -ame ) (18) 


Similar to the anode, the boundary conditions for Eqs. (17) 
and (18) are 


dne fE 
x = 0, Jee = Jiotal, Tx = — Pk c Jiotal; 
dnc og 


x= dg, Jic = Stotal, = pT Jota (19) 


dx 


The transport of O2 can be described by the self-diffusion 
mechanism [20], 


dPo, RTh a ea 
= O2: 02 
dx 4FDĚ P 0 > 
and 
pet 
50, == (21) 


2 = eff eff 
Do, k + Po,- 


The boundary condition for transport of O2 (Eq. (20)) is Po, = 
P, at x=0. The cathode total overpotential can be obtained as 
[13] 


Nuc = boc — Pec) — (Pied) — Pe,ca=0)) (22) 
2.4. The electrolyte 


The dense electrolyte is typically made of yttrial-stabilized 
zirconia (YSZ) and is used to conduct O?~. The overpotential 
of the dense electrolyte can be determined by Ohm’s law, 


Ne = Stotal Re L (23) 


where Re and L are the resistivity (Q m) and thickness (m) of 
the electrolyte, respectively. 


3. Model validation 


In this section, the model as described in the previous section 
is validated by comparison with experimental data. The values 
of the input parameters for subsequent simulation and analyses 
are summarized in Table | [12,13,16,21,22], unless otherwise 
stated. 


Table 1 

Values of input parameters used in the present model 

Parameter Value 

Temperature, T (K) 1073 

Pressure, P (atm) 1.0 

Reference anode exchange current density, 4000.0 [21] 
JEE (Am?) 

Reference cathode exchange current density, 1300.0 [21] 
JE (Am?) 

Anode thickness, da (m) 500 x 1076 

Cathode thickness, de (m) 50 x 1076 

Electrolyte thickness, L (m) 50 x 1076 


Conductivity of electronic conductor of 2x 10° [12:13] 


anode, o¢, (S m7!) 


Conductivity of electronic conductor of 1 x 104 [22] 
cathode, de. (S m7!) 
Conductivity of ionic conductor, © ;,a,0 i,c 3.44 x 


(Sm7!) 10* exp(—10300/T) 
[12,13] 
Electrode tortuosity, € 4.5 [12,13] 
Volume fraction of electronic conductor, ® 0.4 [12,13] 
Electrode porosity, £ 0.3 [12,13] 


Particle size, re, ri (m)* 0.5 x 1076 [16] 


à Note: The values of porosity and particle size are used as typical microstruc- 
ture parameters. Other values for FGE and for comparative analyses are also 
used and specified in the text. 
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Jiang et al. [23] described their laboratory set-up and exper- 
imental test procedures in detail and, therefore, their work was 
selected for model validation. In their experiments, both pure 
Ni and Ni-YSZ cermet anodes with different YSZ loadings 
were tested. The volume percentages of the impregnated YSZ 
in the anode were 17 vol.% and 21 vol.% with YSZ loadings 
of 2.7mgcm~? and 4.0mgcm~?, respectively. In the sample 
preparation, an anode was deposited ona 1-mm YSZ electrolyte. 
The thickness of the anode was about 30 um and its porosity was 
30%. The tests were conducted under a constant pressure of 1 
atmosphere (atm), H? molar fraction of 97% (3% H20), and 
temperature ranging from 973 K to 1173 K. The performance 
of the pure Ni and Ni-YSZ cermet anodes were characterized 
by polarization and electrochemical impedance spectroscopy 
(EIS). Fig. 3(a) shows the dependence of anode overpotential 
on current density at different YSZ loadings. A good agreement 
between the theoretical simulation and the experimental data is 
found. The Ni-YSZ cermet anodes have a lower overpotential 
than the pure Ni anode. It is because the electrochemical reac- 
tions are extended from the EE interface to the inner Ni- YSZ 
electrode, resulting in an enhanced electrode reactivity. Fig. 3(b) 
shows the effect of operating temperature on the overpotential 
of the cermet anode with 2.7 mg cm~? YSZ loading. The anode 
performance increases with increasing temperature due to the 
enhanced ion conduction. Again, the simulated results agree 
reasonably well with the experimental data by Jiang et al. [23]. 

In the second simulation, the J-V characteristics of the SOFC 
were simulated and compared with the experimental data by 
Kim et al. [24]. In Kim et al.’s experiments, the anode, cath- 
ode, and the electrolyte were made of Ni-YSZ, YSZ, and 
Lao.gSro.2MnO3-YSZ, with thickness about 1.2 mm, 30 um, and 
7 um, respectively. Experiments were conducted at a pressure 
of | atm, temperature of 1073 K, and a H? molar fraction of 
97% (3% H20). The measured J-V data were compared with 
the theoretical simulation results as plotted in Fig. 3(c). Again, 
a good agreement between simulation results and experimental 
data was found. Fig. 3(a—c) serve as a thorough validation of the 
present model. 


4. Results and discussion 


In the following analyses, a conventional non-graded SOFC 
and a micro-structurally graded SOFC were studied for their 
electrochemical performance. From the previous studies, SOFC 
electrode particles having radii ranging from 0.012 pum to 
5.0 um have been reported [25-27]. For electrodes made from 
spherical particles, the achievable porosity ranges from 0.3 to 
0.7. For a porosity less than 0.3, the extent of overlap between 
spherical particles increases with decreasing porosity, causing 
the active surface area to decrease considerably with decreasing 
porosity [28-30]. For a porosity higher than 0.7, the electrode 
will suffer from poor particle connectivity, poor percolation, and 
thus poor performance [30,31]. Therefore, for electrodes con- 
sisting of spherical particles, porosity between 0.3 and 0.7 is 
regarded as a practical range. Without losing generality, the most 
favourable anode-supported configuration with thick anode and 
thin cathode/electrolyte is considered in this study. 
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Fig. 3. Comparison between simulation results and experimental data from lit- 
erature for model validation: (a) anode overpotential at different YSZ loadings; 
(b) anode overpotential at different temperatures; and (c) J-V characteristics of 
a planar SOFC. 


4.1. Current density and gas composition distribution in 
electrodes 


In this section, the distributions of electronic/ionic current 
density and gas composition along the electrode depth were 
studied at the typical operating current density of 10,000 A m7. 
Fig. 4(a) shows the distributions of electronic/ionic current den- 
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Fig. 4. Current density and gas composition distribution along electrode depth: (a) electronic and ionic current density distribution in anode; (b) Hz molar fraction 
distribution in anode; (c) electronic and ionic current density distribution in cathode; and (d) O2 molar fraction distribution in cathode. 


sity in porous anode. The electronic current density increases 
from O at the EE interface (x=500 um) to the maximum 
(10,000 A m~’) at the anode surface (x =0), where electrons are 
collected. The ionic current density shows a reverse trend of 
electronic current density. The sharp change of electronic/ionic 
current density at the EE interface indicates that most of the 
electrochemical reactions occur there. It is because the effective 
ionic conductivity is much lower than that of electronic conduc- 
tivity. Fig. 4(b) shows the distribution of H? molar fraction in the 
anode region. The concentration of Hz decreases considerably 
from the anode surface towards the EE interface, indicating the 
importance of gas transport in the thick anode (500 um). The 
electronic/ionic current density distribution in the porous cath- 
ode is shown in Fig. 4(c). The steady changes in electronic/ionic 
current density indicate that electrochemical reactions take place 
throughout the thin cathode (50 um thick). In addition, the thin 
cathode shows a very low mass transfer resistance, as can be 
seen from Fig. 4(d). 


4.2. Effect of electrode microstructure on electrode 
performance 


In order to better understand the effect of micro-structural 
grading, it is important to understand how the micro-structural 


parameters affect the electrode performance. In this section, the 
important electrode porosity and particle size are analysed. 
Fig. 5 shows the effect of particle size on the electrode over- 
potential. For both anode and cathode, optimal particle sizes 
of about 1.0 um and 0.3 wm are found. In general, the smaller 
the particle size, the larger the volumetric reactive surface area 
(indicated by Eqs. (2) and (3)). Therefore, reducing particle 
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Fig. 5. Effect of particle size on electrode overpotential. 
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Fig. 6. Effect of porosity on electrode overpotential. 


size can increase the electrode reactivity and thus reduces the 
electrode activation overpotential. However, when the particle 
size is too small, the reactant concentration at the reactive sites 
will approach to zero and thus the concentration overpotential 
will increase considerably. As the anode in an anode-supported 
SOFC is much thicker than the cathode, the optimal particle size 
of the anode is larger than that of the cathode. 
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Similarly, the effect of electrode porosity is shown in Fig. 6. 
An optimal porosity of about 0.4 is found for the thick anode. 
Reducing porosity can reduce the ohmic overpotential, as more 
solid particles are available to facilitate the transport of elec- 
tronic/ionic charges. In addition, the activation overpotential is 
also reduced with decreasing porosity due to the increase in 
reactive surface area. However, the concentration overpotential 
of the thick anode increases with decreasing porosity since less 
space is available for gas transport. The combined effects of 
porosity on ohmic, activation, and concentration overpotentials 
for minimisation of the overall overpotential result in an optimal 
porosity for the anode. For comparison, the cathode overpoten- 
tial increases monotonically with the increasing porosity. It is 
because the gas transport does not reach its limit in the thin 
cathode of an anode-supported SOFC. 


4.3. SOFC with micro-structurally graded electrodes 


The power density and the electrode overpotentials of SOFC 
with conventional non-graded electrodes and SOFC with poros- 
ity/particle size graded electrodes are compared in Figs. 7-9. As 
anon-graded SOFC can hardly be operated with particle size less 
than 0.3 um due to the large anode concentration overpotential 
(Fig. 5), a linear grading in particle size in the range from 0.3 um 
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to 5.0 um is studied. In addition, a linear grading of porosity in 
the range from 0.3 to 0.7 is also employed. In the present study, 
grading means that the particle size or porosity increases linearly 
from the EE interface to the electrode surface. As the grading 
of electrode does not change the performance of electrolyte, 
the electrolyte ohmic overpotentials are not presented in this 
paper. 

Fig. 7 shows the power density and electrode overpotential 
of SOFCs with large particle size (1.0 wm) at the EE interface. 
The maximum power density of a SOFC with porosity-graded 
electrodes is found to be 23% higher than that of a non-graded 
SOFC (Fig. 7(a)). Grading of electrode porosity can significantly 
reduce the mass transfer resistance, but also reduce the volu- 
metric reactive surface area (Eqs. (2) and (3)), and effective 
electronic/ionic conductivity (Eq. (7)). As the electronic con- 
ductivity is several orders of magnitude higher than the ionic 
conductivity, the ohmic overpotential is primarily caused by 
the resistance to the ionic current, which mainly occurs at a 
thin layer close to the EE interface. In addition, the activation 
overpotentials also mainly take place near the EE interface. For 
a porosity-graded thick electrode, the porosity does not vary 
much in the active thin layer. Thus, the ohmic and activation 
overpotential are insensitive to the porosity grading, especially 
for thick electrodes. As a result, the porosity-graded electrodes 


show smaller overpotential than non-graded electrodes (Fig. 7(c 
and d)). In a previous study conducted by Greene et al. [11], 
the maximum power density of a porosity graded SOFC is only 
about 8.7% higher than that of a non-graded SOFC. There is a 
difference between the present study and Greene et al. [11] due 
to the overestimation of the ohmic overpotential in Greene et 
al.’s study. Their macro-scale model cannot predict the distribu- 
tion of ionic current density and thus they consider the ohmic 
overpotential to occur throughout the electrode, resulting in an 
ohmic overpotential higher than the actual value. The present 
micro-scale model can predict more precisely the performance 
of SOFC with micro-structurally graded electrodes. Compared 
with porosity grading, the particle size grading shows compara- 
ble performance. In addition, the larger the extent of the anode 
particle size grading, the lower the anode overpotential, and thus 
the better the SOFC performance. The grading in cathode can 
also reduce the mass transport resistance and thus shift the limit- 
ing current density to a larger value (Fig. 7(d)). Unlike porosity 
grading, particle size grading does not change the ohmic over- 
potential. However, when the particle size grading is too radical 
for the thin cathode (re: 5.0 um — 10 um), the activation over- 
potential will become more important at low to medium current 
density, as can be seen from Fig. 7(d). Under these circum- 
stances, too much grading of the cathode may result in a poorer 
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SOFC performance (compared with only anode-graded SOFC) 
at medium current density (Fig. 7(b)). 

For particle radius of 0.5 um and 0.3 um at the EE interface, 
the performances of SOFC with or without micro-structurally 
graded electrodes are presented in Figs. 8 and 9, respectively. 
It is found that the smaller the particle size at the EE interface, 
the larger the difference between the micro-structurally graded 
SOFC and the non-graded SOFC (Figs. 8 and 9). For an interfa- 
cial particle size of 0.5 um, the particle size graded SOFC shows 
significantly higher (as much as 70%) power density than the 
non-graded SOFC (Fig. 8(a)). In case of particle radius equal 
to 0.3 um at the EE interface, the maximum power density of 
a SOFC with considerable size grading (r: 5.0 um — 0.3 ppm) 
is 153% higher than that of a non-graded SOFC (Fig. 9(a)). 
For a non-graded SOFC, decreasing particle size results in con- 
siderable increase in the mass transfer resistivity and thus the 
SOFC power density is low due to the small limiting current 
density. As mentioned previously, the concentration overpo- 
tential is more sensitive to the micro-structural grading than 
the ohmic and activation overpotentials. Since the gas trans- 
port becomes more important for a SOFC with small particle 
size than one with large particle size, the performance gain 
is more pronounced for micro-structurally graded SOFC due 
to the considerably decreased mass transport resistance (Figs. 


8(b and c), and 9(c and d). In addition, it can be seen that a 
large extent of grading of the thin cathode is needed, as illus- 
trated in Fig. 9(b) that a SOFC with both electrodes graded has 
a higher power density than a SOFC merely with the anode 
graded. 

Compared with porosity grading, particle size grading seems 
to be more effective to enhance the SOFC performance for small 
EE interfacial particle size. One reason is that the porosity of 
the electrodes can only change in a small range (0.3-0.7) due 
to the limitation of percolation [10,29]. Another reason is that 
the ohmic overpotential will increase due to grading in poros- 
ity. On the contrary, particle size grading does not alter the 
ohmic overpotential. More importantly, the particle size can 
be changed in a wider range. As a result, grading in particle 
size is superior to porosity grading since the current research 
trend is to make particle size smaller and smaller at the EE 
interface to enhance the electrode reactivity at an intermediate 
temperature [8,9,25—27]. In fuel cells and other fields, micro- 
structurally graded electrodes can be fabricated by a number of 
ways, such as combustion chemical vapour deposition (CCVD) 
[8] and electrophoretic deposition [32]. Therefore, the concept of 
micro-structurally graded SOFC is feasible and has great poten- 
tial to improve the SOFC performance by simply modifying the 
electrode microstructures. 
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5. Conclusion 


A mathematical model has been developed to investigate the 
effect of micro-structural grading on the SOFC performance. 
The model was validated by comparing the simulation results 
with experimental data available in the literature. The model 
is able to analyse the coupled phenomena of electrochemical 
reactions and mass transfer in the SOFC electrodes. 

Itis found that both porosity grading and particle size grading 
are effective to enhance SOFC performance. With a decrease 
in particle size at the EE interface, the performance gain by 
particle size grading method becomes more pronounced. It has 
been demonstrated that the maximum power density of a parti- 
cle size graded SOFC could be 153% higher than a conventional 
non-graded SOFC with small particle size, indicating the great 
promise of micro-structural grading in SOFC. The reason is that 
the concentration overpotential can be considerably reduced by 
proper grading of electrode microstructure while the ohmic and 
activation overpotentials are changed only slightly with micro- 
structural grading. It is found that the particle size grading 
becomes more advantageous than porosity grading when the 
particle size at the EE interface is small. It is also found that 
the effect of grading is generally beneficial for thick electrodes. 
However, when the EE interfacial particle size is large, too much 
grading in particle size of the thin cathode should be avoided for 
SOFC operating at a low to medium current density. 

This paper reveals that micro-structural grading is very 
promising to enhance the SOFC power density and worthy 
of further investigation. Aside from their electric performance, 
other properties associated with micro-structural grading, such 
as mechanical strength and catalyst sintering, should be further 
investigated. 
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